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The state-of-the-art in vehicle design decouples flight feasible trajectory generation from 
the optimization process of an entry spacecraft shape. The disadvantage to this decoupled 
process is seen when a particular aeroshell does not meet in-flight requirements when 
integrated into Guidance, Navigation, and Control simulations. It is postulated that the 
integration of a guidance algorithm into the design process will provide a real-time, rapid 
trajectory generation technique to enhance the robustness of vehicle design solutions. The 
potential benefit of this integration is a reduction in design cycles (possible cost savings) 
and increased accuracy in the aerothermal environment (possible mass savings). This work 
examines two aspects: 1) the performance of a reference tracking guidance algorithm for five 
different geometries with the same reference trajectory and 2) the potential of mass savings 
from improved aerothermal predictions. An Apollo Derived Guidance (ADG) algorithm 
is used in this study. The baseline geometry and five test case geometries were flown 
using the same baseline trajectory. The guided trajectory results are compared to separate 
trajectories determined in a vehicle optimization study conducted for NASA’s Mars Entry, 
Descent, and Landing System Analysis. This study revealed several aspects regarding the 
potential gains and required developments for integrating a guidance algorithm into the 
vehicle optimization environment. First, the generation of flight feasible trajectories is 
only as good as the robustness of the guidance algorithm. The set of dispersed geometries 
modelled aerodynamic dispersions that ranged from +/-1% to + /-17% and a single extreme 
case was modelled where the aerodynamics were approximately 80% less than the baseline 
geometry. The ADG, as expected, was able to guide the vehicle into the aeroshell separation 
box at the target location for dispersions up to 17%, but failed for the 80% dispersion cases. 
Finally, the results revealed that including flight feasible trajectories for a set of dispersed 
geometries has the potential to save mass up to 430 kg. 

I. Introduction 

Advances in optimization algorithms and computational capabilities have lead to the development of 
Multi-Disciplinary Design, Analysis, and Optimization (MDAO) environments for entry spacecraft design. 
MDAO environments integrate models from multiple fields of interest to generate candidate vehicle con- 
cepts that meet given mission criteria. The objective functions in entry vehicle optimization can be a 
function of trajectory dependent variables, for example, minimizing peak heat rate where heat rate is a 
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function of the altitude-velocity condition. The MDAO environment may include parametric vehicle shapes, 
trajectory optimization, structures models, and thermal protection system (TPS) sizing models. The ae- 
ordynamic/aerothermodynamic databases for the ’optimal’ candidate geometry are provided to Guidance, 
Navigation, and Control (GNC) engineers who develop relevant guidance and control schemes to achieve 
trajectory specific mission criteria. GNC engineers complete dispersion analyses to test the performance of 
GNC algorithms in meeting targeting requirements. The dispersions include, but are not limited to, uncer- 
tainties in navigational position, small changes in entry interface conditions, and uncertainty in aerodynamic 
quantities. The dispersion analysis generates a large database of heating and loading profiles corresponding 
to each dispersed trajectory. This database of trajectories is passed back to the MDAO tool to improve the 
candidate vehicle shape and TPS configuration. 

The trajectory model in most MDAO environments does not necessarily determine a new trajectory when 
the geometry changes. In addition, the trajectory model does not always guarantee that a flight feasible 
trajectory is used to determine trajectory dependent parameters, such as heat rate. Herein, a flight feasible 
trajectory is defined as a profile that obeys the control limits and meets the targeting requirements. The use 
of trajectories that are not updated or not flight feasible may lead to an over-designed TPS configuration 
and may not reveal critical design errors that are eventually found after testing a physical prototype. 

This research investigates the performance of an Apollo-Derived Guidance 1 algorithm in the presence of 
dispersed aerodynamics due entirely to changing the entry vehicle geometry, which is what occurs in a shape 
optimization study. The use of a guidance algorithm as a trajectory model provides the potential to rapidly 
update the trajectory for a new geometry and ensure that a flight feasible trajectory is generated. 

I. A. Background 

There are multiple shape optimization studies that have been completed, where some may or may not include 
trajectory models or trajectory dependent objective functions. In Theisinger et al. a set of analytical 
geometric functions were implemented into an optimizer to study multiple shapes that minimize the ballistic 
coefficient while maintaining a specific L/D. However, no trajectory dependent variables were considered 
in the shape optimization, e.g. heat rate and/or heat load. There are some shape optimization studies 
that use aerothermodynamic analyses to optimize shapes. In Johnson et al." , the shape of hypersonic 
blunt-body heat shields were optimized by minimizing total heat rate at the stagnation point. The vehicle 
speed and altitude are based on a single representative trajectory, thus flight feasible trajectories were 
not determined as the shape changed in the optimization. Another MDAO tool, called HAVOC 4 , was 
used to propose a conceptual design for a TPS material called Ultra High Temperature Ceramics (UHTC). 
HAVOC contained a 3DOF trajectory model consisting of two methods to compute the trajectories. The 
first is an energy state approximation method coupled with a Mach number vs. altitude scheduling. The 
second method solves the equations of motion with 1st order integration subject to a user specified set 
of path constraints. This tool provided engineers with the ability to complete a number of design cycles 
such that an optimum design was determined. However, the trajectory analysis lacked flight-feasible bank 
profiles and thus does not account for scenarios which may adversely affect the vehicle shape and TPS 
design. In Garcia et al. 5 an MDAO tool, called Co-optimization Bluntbody Re-Entry Analysis (COBRA), 
was used to determine optimal aeroshell shapes for Mars atmospheric entry. COBRA includes parametric 
geometry models, aerodynamic models, aerothermodynamics models, trajectory optimization, TPS analysis 
tools, structures analysis, and mass estimation models. The objective functions were to minimize peak 
heat rate and maximize CD A to minimize the ballistic coefficient for a mid-L/D shape. While this shape 
optimization includes almost all disciplines including trajectory optimization, the trajectory analysis does 
not necessarily determine flight feasible trajectories. In Steinfeldt et al. 6 , a vehicle architecture landing 
assessment is completed with the integration of, among other models, a trajectory model. It is specifically 
stated in this paper that a close-loop guidance algorithm was not be included in the analysis because it would 
constrain the predicted performance of the system architecture. However, it can be argued that modelling 
closed loop guidance will provide the engineer with a realistic study of the performance that is needed to 
achieve the mission objectives. In Loomis et al. and Saunders et al. 8 , a TPS sizing tool was developed 
that integrates aerodynamics, aerothermodynamics, and trajectory analysis. In Loomis' the size of a TPS 
aeroshell for the Mars Smart Lander was baselined. The authors use trajectory optimization and TPS sizing 
to achieve this. A constraint on the distance to the target landing site is included, with the potential to 
include this value in the objective function. However, the study conducted focuses on maximizing cross 
range instead of achieving a precision landing. In addition, there does not appear to be any constraints 
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on control rate or control acceleration and does not indicate how long this simulation takes to converge 
to a solution. Another tool for optimizing vehicle shapes is found in two papers authored by Kinney et 
al.' ' ’ . The Configuration Based Aerodynamics (CBAERO) tool was used to optimize the HL-20 vehicle 
shape by maximizing L/D while holding volume constant and constraining the pitching moment to zero. In 
this study, no trajectory dependent parameters were used as constraints or objective functions. In Johnson 
et al. 11 , a population-based differential evolutionary algorithm was used to determine the optimal vehicle 
design for the multi-objective function to minimize heat load and maximize cross range. The class of vehicles 
under consideration ranged from an L/D of 0.3 to 1. The authors selected entry conditions similar to the 
Orion spacecraft and used POST to model the trajectory. The bank profile was assumed to have a step 
function behavior. An initial guess of the bank profile was input and then bank angle adjustments were 
made through trial and error, until a trajectory was generated to meet the constraints. The trajectories 
were not necessarily flight feasible trajectories since targeting and bank limits (rate/acceleration) were not 
included. Clearly many MDAO studies lack the verification of flight feasible trajectories or any update 
once the trajectory changes. This gap in trajectory modelling for entry vehicle shape optimization is the 
motivation for studying guidance performance in the presence of this type of aerodynamic dispersion. 

II. Method 

As a pre-cursor to full integration, this study aims to determine the potential advantages to integrating 
a guidance algorithm into an entry vehicle MDAO environment. It is not within the scope of this study to 
develop a guidance algorithm or complete a full entry vehicle analysis. The general test case for this study 
is a Mars direct entry of a rigid mid-L/D aeroshell to land a 40 metric ton payload. 

II. A. Guidance and Trajectory Modeling 

In order to understand the limitations and/or potential of integrating current guidance algorithms into an 
MDAO environment, the Apollo Derived Guidance (ADG) algorithm was selected for this study. The ADG 
was used for Apollo’s atmospheric entry and is widely used by the entry spacecraft GNC community. The 
ADG is a reference tracking algorithm with a terminal point controller derived using the calculus of variations 
for linear perturbations. The reference trajectory consists of drag acceleration, altitude rate, range to go, 
vertical L/D, and controller gains as a function of relative velocity. The application of the ADG in this work 
will provide an initial understanding of guidance performance in the presence of aerodynamic dispersions 
that are due entirely to a change in geometry. Although trajectory design is outside the scope of this study, 



Figure 1. A) Reference drag acceleration, B) altitude rate, 



C) range to go, and D) vertical L/D versus relative velocity. 


selecting an appropriate reference trajectory for implementation into the ADG is critical. The reference 
trajectory must ensure that the vehicle reaches its target and that the control profile does not violate control 
system constraints. In addition, the trajectory must ensure that aerodynamic loading and heating constraints 
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are not violated. A nominal trajectory for a baseline mid-L/D vehicle shape was designed by NASA Langley 
Research Center that meets these criteria. This trajectory and corresponding geometry can be found in 
NASA’s comprehensive Mars Entry, Descent, and Landing Systems Analysis. 1 This baseline geometry and 
trajectory was implemented into the ADG for this work. This baseline is well suited for this study since the 
results will be compared with an entry vehicle MDAO study that used the same baseline case. The reference 
trajectory variables for the ADG are plotted in Figure 1. 

The guidance algorithm is run in the Simulation of Rocket Trajectories (SORT) program. 1 MarsGRAM 
2005 is used to model the Mars atmosphere, the gravity model accounts for J2, J3, J4 -harmonics, and 
the aerodynamic and aerothermodynamic databases are acquired from CBAERO. The databases are a 
function of Mach number, dynamic pressure, and angle of attack, where angle of attack is fixed for each 
vehicle. The termination criteria (Table 1) correspond to separation of the aeroshcll for the descent phase. 
In addition, the aerodynamic loading and heating constraints are monitored for each trajectory according to 
specified limits (Table 1). 


Separation 

h[km\ 

q[N/m 2 ] 

Mach// 

Lower 

6.4 

1077.7 

2.4 

Upper 

8.9 

1598.3 

3.3 

Maximum 

q[N/m 2 ] 

Loading[g’ s\ 

q'\W/cm 2 ] 


15000 

4.0 

1000.0 


Table 1. Trajectory parameters and termination criteria. ' 


II. B. Vehicle Shapes and Aerodynamic Dispersions 

The total entry mass for the baseline geometry is 109.59 metric tons corresponding to a ballistic coefficient 
of 401.7 kg/m 2 . Figure 2 illustrates the geometry with speed breaks and flaps, however, these control surface 
deflections were set to zero to generate the baseline aerodyamic and aerothermodynamic databases. This 
geometry will be referred to as the baseline or ’sled’ throughout this paper. From this point forward, dispersed 


Side View 



Figure 2. Baseline Sled Geometry 


geometries will refer to 5 varied geometries generated from a comprehensive MDAO study of rigid, mid-L/D 
aeroshells. The target L/D for all geometries is 0.5, where the change in geometry results in a change in the lift 
and drag coefficients but not the resulting ratio. In addition, the reference area and reference length for four 
of the five dispersed geometries is 78.54 square meters and 10 meters, respectively. In Garcia ’, the aeroshcll 
design space was explored by wrapping a Multi-Objective Genetic Optimization around an integrated set of 
relevant models. The Co-Optimization of Blunt-body Re-entry Analysis (COBRA) technique parametrically 
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changes the aeroshcll shape and uses a single nominal trajectory to define the aerodynamics loading and 
heating characteristics. Essentially, a single peak heating ’altitude-velocity’ condition corresponding to a 
single reference geometry, is used to calculate the heat rate for each candidate geometry. The point of 
departure (POD) geometry, which is a spherical nose attached to a cylinder (Figure 3A), is the starting 
point for the parametric changes. The objective functions were to maximize CdA and minimize peak 
heating. The constraints included the launch shroud envelope and static stability for pitch and yaw. The 
resulting pareto front from the genetic optimization is a scatter plot of the total peak heat rate and CdA 
for each geometry in the design space. In the COBRA study, three candidate vehicles were selected such 
that the trade off between heat rate and CdA were considered. In Garcia ', these three candidate geometries 
were then run through the COBRA process again, but with some modifications. First, a pre-processed 
trajectory is optimized to match the altitude-velocity profile of the baseline trajectory. This trajectory does 
constrain the aerodynamic loading and heating, however, there are no constraints on the controls. Also, 
short of targeting a particular termination longitude and latitude for the direct entry portion of flight, the 
trajectory was constrained to ensure that the downrange was larger than 1220 km. The trajectory analysis 
for the COBRA work included aerocapture and direct entry. The study in this paper focuses on the direct 
entry for the dispersed geometries. Second, a mass closure model iterates on the total entry mass until the 
subsystem masses (structure, TPS, etc.) converges. At this point a comparison of the aerodynamic loading, 
heating characteristics and total mass is completed to determine the ’’optimal” vehicle design. Four of the 
five dispersed geometries used in this study include the POD and the three candidate geometries selected 
from the COBRA study. It should be noted that the aerodynamic dispersion is a deviation from the Cl and 




Aspect Ratio 
A through D: 10x30 meters 
E: 10x24 meters 


Figure 3. Dispersed geometries from COBRA study. 


Cd used to generate the baseline trajectory. The first dispersed geometry, COBRA POD, has a ballistic 
coefficient of 492.1 kg/m 2 and a wetted area of 1021 meters squared. The dispersion on the drag and lift 
coefficients is -7.5% and -13.3%, respectively, from the baseline case. The next dispersed geometry is COBRA 
8459, which has the lowest CdA and lowest heat rate value on the pareto front. The ballistic coefficient is 
406.04 kg/m 2 and the wetted area is 877 meters squared. The dispersion on the drag and lift coefficients is 
+4.2% and -1.4% respectively. The third dispersed geometry is COBRA 14297, which has a CdA and heat 
rate that is mid-range on the pareto front. The ballistic coefficient is 365.7 kg/m 2 and the wetted area is 
882 meters squared. The dispersion on the drag and lift coefficients is +15.8% and +9% respectively. The 
fourth dispersed geometry is COBRA 14888, which has the highest CdA and highest heat rate on the pareto 
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front. The ballistic coefficient is 368.03 kg/m 2 and the wetted area is 937 meters squared. The dispersion 
on the drag and lift coefficients is +17.7% and +11.15% respectively. In the COBRA study it was found 
that COBRA 14297 had the optimal heating and weight characteristics. To test the ADG performance for 
an extreme case, the length of this vehicle was shortened and the aerodynamics were scaled to reflect the 
change in shape. This case is referred to as Scaled Aspect Ratio (AR) throughout this paper and defines 
the fifth dispersed geometry for this study. The ballistic coefficient is 922.4 kg/m 2 and the wetted area is 
364.21 meters squared. The dispersion on the drag and lift coefficients is -70% and -85% respectively. All 
geometries can be found in Figure 3. 


III. Results 

Key questions must be addressed in analyzing the results of this study. 

• Did the vehicle meet the separation box defined by Mach Number and Dynamic Pressure? 

From the perspective of entry vehicle analysis, it is important that the vehicle reach the separation 
box. The box is defined by the Mach number and dynamic pressure limits found in Table 1. If the 
control profile does not take the vehicle within the separation box then the control profile is invalid. In 
Table 2 the separation condition for each vehicle from the COBRA study and the ADG are compared. 
The COBRA 8459, 14297, and 14888 are sufficiently guided into the separation box by the ADG. The 
ADG is able to guided the COBRA POD close to the separation box, where the dynamic pressure is 
approximately 28 N/m 2 larger than the upper limit. The ADG was not able to guide the Scaled AR 
into the separation box as the dispersions on the lift and drag coefficients were too large. 


Separation 

h[km] 

q[N/m 2 ] 

Mach // 

Nominal 

6.41251 

1232.572 

2.67 



COBRA Study 

Guided Study 


h[km\ 

q[N/m 2 ] 

Mach// 

h[km] 

q[N/m 2 } 

Mach// 

POD 

5.8525 

1432.69 

2.81 

6.41471 

1626.142 

3.06 

8459 

6.2630 

1259.06 

2.67 

6.41810 

1320.308 

2.76 

14297 

6.1497 

1272.87 

2.67 

6.41139 

1250.032 

2.69 

14888 

6.0297 

1287.27 

2.67 

6.41139 

1250.908 

2.70 

Scaled AR 

5.9920 

1269.39 

2.65 

6.41813 

28665.294 

12.85 


Table 2. Separation condition comparison between the COBRA study and Guided results. 


As with any trajectory optimization, it is possible to always find a control profile that takes the vehicle 
to a given separation condition, which is the case for the COBRA study trajectories. In Table 2 all 
the COBRA trajectories are shown to meet the separation box. However, there are no constraints that 
identify whether any control limits were violated or if the vehicle is within a reasonable radius of the 
target. The worst case scenario is that the optimized trajectory is not a flight feasible trajectory and 
thus does not model the heating environments accurately. 
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Bank Command |dcg] Bank Command [deg| Bank Command [deg] Bank Command [dcg| Bank Command |dcg[ 
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Crossrange From Target |nmi| 


In Figure 4, the bank profile for each dispersed geometry from the COBRA study and the ADG are 
plotted with respect to the nominal trajectory. It is noted in Mendeck and Craig 15 that the ADG 
includes a limiter on the commanded vertical L/D such that bank cannot exceed +/- 90 degrees. The 
purpose of this limiter is to ensure that the g-loading does not become excessive and violate loading 
constraints. This condition is especially critical for Mars entry because a full lift-up orientation (0°) 
for the duration of entry will guide the vehicle to a separation state well beyond the target and a 
full lift-down orientation (180°) will not allow adequate slowing of the vehicle, leading to catastrophic 
surface impact significantly short of the target. In Figure 4 the guided bank profiles (solid blue lines) 
for all dispersed trajectories are within the +/- 90 degree limits. The COBRA Study bank profiles 
(blue dashed lines) for the POD, 8459, and 14888 do not fly within the required bank limits. Thus, it 
is possible that the heating environments modelled do not reflect the actual in-flight behavior. Using 
the COBRA trajectory optimization process, a bank profile for the Scaled AR was found that took the 
Scaled AR to the separation box, but did not remain within the limits of +/- 90 degrees. The ADG, 
however, flew a bank solution that maintained the control limits, but could not take the vehicle to the 
target. This results lead to the next question to consider. 

• Does the prescribed bank profile guide the vehicle to within a reasonable distance of the target? 

The goal of integrating a guidance algorithm in a shape optimization study is to model a realistic set 
of controls. This implies that the vehicle flics a realistic profile that takes the vehicle to its intended 
destination and does not exceed control limits. The target in this case is the separation condition for 
the descent stage, which is approximately 22.9 km from the target. This separation condition will 
become the ’zero’ range for all trajectories and any deviation from this point indicates the amount 
of range error present from the separation condition. Figure 5 plots the ’separation’ footprint for all 
trajectories considered. It should be noted that the cross range for the COBRA study trajectories 
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Figure 5. Separation footprint for all geometries and corresponding trajectories. 


was set to zero because it was not calculated with respect to the target condition and should not be 
mistaken for a targeted condition. The guided trajectories show moderate performance. The COBRA 
14297 dispersed geometry is guided by the ADG to within 2 km of the location of separation. The 
COBRA 8549 and 14888 vehicles were guided to within 5km of the target and the COBRA POD is 
just outside of 30 km from the target. The case with the largest violation is the Scaled AR, which 
is far outside the target. Note that all of the trajectories from the COBRA study have considerable 
downrange miss distance, between 222 km and 257 km. At this point, one must determine whether 
reaching the separation box far from the target location versus targeting both the separation box and 
target location dramatically changes the heating & loading environment predictions. 

• What are the differences in the aerodynamic loading and heating characteristics between each trajectory 
model? 

At this point the peak conditions for each dispersed geometry and corresponding trajectory are com- 
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pared. However, the Scaled AR is not considered from this point on as it is clearly not a valid case. 


Maximum 

q [N/m 2 ] 

g-load [Gs] 

q [W/cm 2 ] 

Heatload [J/cm 2 ] 

Nominal 

11642.27 

3.25 

125 

12851.97 



COBRA Study 


q [N/m 2 ] 

g-load [Gs] 

q [W/cm 2 ] 

Max. Heatload [J/cm 2 ] 

POD 

10896 

2.59 

89.495 

9445.8 

8/59 

10619 

2.89 

84.727 

11496.4 

1/297 

9555.8 

2.92 

81.7354 

11609.5 

1/888 

9456.1 

2.92 

82.524 

11669.1 


Guided Study 


q [N/m 2 ] 

g-load [Gs] 

q [W/cm 2 ] 

Max. Heatload [J/cm 2 ] 

POD 

13347.20 

3.19 

97.924 

9175.6 

8/59 

11117.93 

3.26 

65.986 

10797.8 

1/297 

10041.16 

3.26 

69.025 

11114.9 

1/888 

10066.52 

3.25 

70.382 

11419.1 


Table 3. Maximum loading and heating comparison between the COBRA study and ADG results. 


In Table 3, the dynamic pressure and g-load for all four dispersed geometries, are higher for the ADG 
trajectories versus the COBRA trajectories. The peak heat rate for the COBRA 8459, 14297, and 
14888 is larger for the COBRA study trajectories as compared to the ADG trajectories. But the 
COBRA POD peak heat rate for the ADG trajectory is larger than the COBRA study trajectory. The 
COBRA POD does not follow the same heat rate trend because this geometry and trajectory were 
calibrated and optimized as the new baseline for the COBRA vehicle analysis, where the basis for the 
new baseline is the Langley trajectory. 

The overall trends that are observed can be attributed to certain aspects of the COBRA trajectory 
optimization process. First, the COBRA trajectories were optimized with an active dynamic pressure 
and g-load constraint. This contributed to control solutions that flew trajectories with excess range 
and excess time in-flight. In Figure 6 the extended time of the COBRA trajectories are 75-100 seconds 
longer than the nominal trajectory. Since the ADG trajectories must reach the same separation box 
and target a specific location, the in-flight time is reduced but the peak dynamic pressure & peak g-load 
are slightly larger. In addition, the entry interface conditions were allowed to moderately change for 
the COBRA trajectory optimization. Specifically, the entry flight path angle (FPA) was made slightly 
steeper by approximately 1.1 degree and the entry altitude was raised by 4 km. Since entry trajectories 
are very sensitive to changes in FPA, the over all result was an increase in the peak heat rate for the 
COBRA trajectories. The final variable for comparison is the integrated heatload, which is a primary 
driver to the size and mass of the thermal protection system (TPS). Heatload is the total heat rate 
integrated over time and is directly proportional to the thickness of the TPS. In Table 3, the heat 
load is shown to be larger for the COBRA trajectories versus the ADG trajectories. This is due to 
the extended trajectory times and higher peak heat rates for the COBRA trajectories. To illustrate 
the effective change in mass from lower predicted heatloads, consider a TPS for this vehicle that is an 
ablator. In Laub et al. L6 a relationship between the TPS mass fraction and the integrated heatload 
was extrapolated from past entry missions that utilized an ablative TPS. The relation is as follows: 

TPS MF = 0.091(Q) 0 ' 51575 (1) 

where MF is the % mass fraction and Q is the total integrated heatload in J/cm 2 . The TPS MF for 
COBRA 8459, 14297, and 14888 is between 0.13 and 0.36% less for the ADG trajectories as compared 
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to the COBRA trajectories. Given a 40 mT landed payload, the total entry mass to land that payload 
may be as high as 120 mT. Thus, the potential savings from using ADG trajectories to size the TPS 
is between 150 and 431 kg. 

IV. Conclusion and Future Work 

The results revealed that including flight feasible trajectories for a given dispersed geometry has the 
potential to save mass up to 430 kg. In addition, this study also revealed that the generation of flight 
feasible trajectories is only as good as the performance of the guidance algorithm. It is well known that the 
ADG is only valid for short range, low to mid-L/D entry spacecraft performing a direct entry. In addition, the 
ADG will perform well as long as the dispersed conditions are within the linear perturbations used to derive 
the gains for the guidance law. In this study, the dispersed geometries modelled aerodynamic dispersions 
that ranged from +/-1% to +/- 17 % and a single extreme case was modelled where the aerodynamics were 
approximately 80% less than the baseline geometry. This algorithm, as expected, was able to guided the 
vehicle into the separation box at the target location for dispersions up to 17%, but failed for the 80% 
dispersion cases. This pre-cursor study has provided a pathway for future work that will be required to 
integrate guidance into the vehicle optimization process. First, a study should be completed to determine 
how guidance performance varies for different baseline trajectories/geometries for a given vehicle and mission. 
Next, all of the geometries that generated the pareto front in the COBRA analysis should be implemented 
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into the ADG so a full statistical analysis is completed to define the ADG performance limitations. Finally, 
an investigation of the procedures required to generate reference trajectories, for each set of geometries, that 
are outside of aerodynamic dispersion limits for the ADG should be completed. 
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